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ABSTRACT  DGPS  system  principles  are  introduced.  Analyses  and 
discussions  are  made  of  the  major  residual  errors  in  carrier  wave 
phase  DGPS  positioning  and  influences  under  low  dynamic  states.  A 
type  of  algorithm  is  put  forward  which — in  carrier  wave  phase  DGPS 
receivers — does  not  require  solving  for  integral  cycle  ambiguity 
but  calculates  positioning  rapidly. 

SUBJECT  TERMS  tCarrier  phase  DGPS  +Fast  positioning 
+Measuring  sequence  tWindow  type  computation 


1  INTRODUCTION 


Following  along  with  the  establishment  of  global  satellite 
systems,  the  GPS  satellite  positioning  system  has  already  far,  far 
surpassed  the  original  intentions  of  the  U.S.  Defense  Department. 
Not  only  has  it  set  off  an  upsurge  of  research  applications  in 
international  armament  technology,  it  has,  moreover,  achieved  wide 
spread  applications  in  a  great  many  civilian  lines  of  work  as  well. 

As  far  as  difference  GPS  (DGPS)  technology  is  concerned — 
particularly,  carrier  wave  phase  DGPS  technology — during  medium  and 
low  dynamics,  high  precision  measurement  positioning,  due  to 
limitations  imposed  by  U.S.  government  secrecy  associated  with  P 
code  (Y  code) ,  it  possesses,  moreover,  relatively  higher  precison. 
Therefore,  it  has  even  broader  prospects  in  various  national 
military  and  civilian  lines  of  work. 

2  DGPS  DETERMINATION  PRINCIPLES 

DGPS  technology  is  to  make  use  of  correlations  associated  with 
GPS  positioning  errors,  carrying  out  mutual  reductions  and 
elminating  common  conjugate  error  terms  in  order  to  increase 
precision. 

^  DGPS  systems  include  a  reference  station  and  several  user 
stations.  Through  taking  GPS  receiver  positioning  results  and  the 
coordinates  of  locations  that  are  already  precisely  known, 
reference  stations  do  comparisions  to  obtain  real  time  positioning 
correction  quantities.  In  conjunction  with  this,  they  take  these 
correction  quantities  and  broadcast  them  out  using  radio  waves. 
Then,  a  method  is  used  where,  through  using  received  positioning 
correction  quantities  broadcast  by  reference  stations,  user 
stations  do  real  time  corrections  on  user  GPS  receiver  positioning 
results  in  order  to  increase  GPS  positioning  precision.  This  is 
nothing  other  than  the  basic  principle  associated  with  GPS. 

General  correction  methods  include  direct  corrections  and  indirect 
corrections . 

(1)  Direct  Correction  Methods:  Reference  station  X,  Y,  Z 
coordinate  correction  quantities  which  are  calculated  out  are 
directly  added  into  GPS  receiver  positioning  solutions.  The 
drawback  of  this  method  is  that,  following  along  with  increases  in 
the  distance  from  user  GPS  receivers  to  reference  stations, 
improvements  in  positioning  accuracy  will  clearly  diminish. 

(2)  Indirect  Correction  Methods:  First,  calculate  out  the 
pseudo  range  correction  amount  associated  with  each  GPS  satellite 
tracked  by  reference  stations.  Then,  transmit  to  users.  User  GPS 
receivers  carry  out  calibration  with  respect  to  pseudo  ranges 
before  doing  positioning  calculations.  An  introduction  is  made 
below  of  the  positioning  correction  principles  associated  with  this 
method. 

Assume  that  the  pseudo  range  from  ground  reference  station  R 
to  the  ith  GPS  satellite,  that  is,  the  pseudo  range  obtained  for 
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the  ith  reception  channel,  is 


/22 


p'=  p\+  (Cdi  -  dT'’)+  dp‘  +  dp  . ' ’  +  d^  ' '+  dM  +  K 
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— true  distance  from  reference 
station  R  to  the  ith  satellite 
— clock  deviation  of  reference 
station  R  relative  to  GPS  time 
system 

— speed  of  light,  m/s 

— ith  satellite  clock  deviation 
relative  to  GPS  time  system 
— range  error  caused  by  GPS  ephereris 
error  (includes  artificial  SA 
errors) 

— range  error  caused  by  troposphere 
time  delay 

— range  error  caused  by  ionosphere 
time  delay 

— range  error  caused  by  receiver 
multiple  path  effects 
— range  error  caused  by  receiver 
measurement  noise 


From  already  known  values  associated  with  reference  station 
three  dimensional  coordinates  and  GPS  ephereris,  it  is  possible  to 


calculate  P„''  .  Moreover,  p\\  can  be  measured  by 

'  I 

reference  station  receivers.  From  equation  (1) ,  one  obtains  the 
psuedo  range  correction  quantity  Wr  i 


I  I  I  - 

^Pr^Pr-Pr 

=  -  C(dT^-  dT^)-  dp;-  dp„^;-  dpj^- dM^- 


(2) 


With  regard  to  pseudo  ranges  measured  by  user  receivers,  from 
equation  (1),  one  has 


p\=  pL=  -  dT;)+  dp;+  dp  dp  '  +  dM^.+  V 


(3) 
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If  use  is  made  of 


to  carry  out  user  range  measurement 
corrections,  then,  from  equations  (2)  and  (3) ,  one  has 

Ap;+  p[=  p'+  C{dx^- dT,)+  (dM,-  (d^,)+  (K,-  V^)  I 

=  (dp^-  dp')+  (dp,.,;  -  dp  J)+  (dp*„^-  dp„,^;)  I 

When  distances  between  GPS  receivers  and  reference  stations 
are  relatively  close  (100km) 


dp'^=  dp^i  ^Pkion~  ^Prien  '  ^P  rion  '^P  rirep 


(4) 

If  reference  stations  and  users  simultaneously  observe  the 
same  4  satellites,  then,  one  has  4  simultaneous  equations  of  the 
type  of  equation  (4) .  From  these,  it  is  possible  to  solve  for  the 
unknowns  Xk,  Yk,  Zk  and  Adri  ,  Adri  to  include  various  residual 
error  terms  obtained  by  the  ith  reception  channel  associated  with 
the  same  observation  epoch 

Adr'=  C  (dij^-  dT,)+  (dMj-  dJl/,)+  V 

3  CARRIER  WAVE  PHASE  DIFFERENCE  GPS  FAST  POSITIONING  CALCULATIONS 
With  regard  to  carrier  wave  phase  observation  quantities  /23 

p[=  ^  (^L+ K)+it>['  (5) 

In  the  equation:  n‘  — is  initial  integral  cycle 

ambiguity 

N'^  — is  integral  cycle  change  number 

beginning  from  observation  epoch 
start  instant  to 

<I)'  — is  observation  phase  decimal 

*  portion 


Therefore , 


Ap;+  p;=  p;+  c(dT^-  dT^)+  (dM,-  dM^)+  (f^-  f,) 
=  [A^'-  Ar*)^+  {Y‘  r*)^+  (Z'-  2*)-]'''^+  Adr' 
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X 


— is  carrier  wave  wave  length. 

Taking  equations  (2)  and  (5) ,  respectively,  and  using 
reference  station  and  user  device  expressions  to  substitute  into 
equation  (4) ,  one  obtains  the  carrier  wave  phase  difference  GPS 
observation  ecpiation 

Pri+  <)+  ^(K+  <)+  K-  ^‘r  i 

=  [(jr'-  (7‘-  (z'-  z*)-]'^"+  Adr''  ‘ 


(6) 

During  general  carrier  wave  phase  measurements,  the  speed 
associated  with  initial  integral  richness  is  the  key  to  increasing 
measurement  speeds.  Classical  measurement  methods  are  to  opt  for 
the  use  of  stationary  or  medium  and  low  dynamics  measurements  made 
for  long  time  periods  on  a  base  line.  Making  use  of  a  batch  of 
observation  data  obtained — at  the  same  time  as  calculating  out 
common  unknown  parameters — calculations  are  made  for  many  estimate 
values  associated  with  initial  integral  cycle  precision  values. 
From  error  ranges  associated  with  the  calculation  values  in 
question  determination  is  made  of  the  only  whole  number  to  act  as 
the  initial  integral  cycle  precision  value.  As  far  as  new  methods 
are  concerned — for  example,  the  Leica  Company's  GPS  measurement 
technology  which  does  not  have  an  initial  common  dynamic  state — 
there  is  still  a  need  to  make  measurement  sequences  of  over  200s. 
At  the  same  time  as  location  coordinates  are  calculated,  option  is 
made  for  the  use  of  fast  approximation  techniques  and  common 
initial  integral  cycle  precision  values  with  uninterrupted 
accuracy,  thereby  obtaining  the  only  exact  initial  integral  cycle 
precision  value  to  act  as  initial  integral  cycle  precision  value 
for  later  measurements  and  accurate  positioning.  In  conjunction 
with  this,  the  methods  in  question  require  that  the  distances 
between  DGPS  reference  stations  and  user  receivers  be  within  10km. 
Among  the  important  causes  for  this  is  that,  due  to  the  fact  that 
when  reference  stations  and  user  receivers  are  separated  by  long 
distance  differences,  the  influences  of  residual  errors  lead  to  an 
inability — within  relatively  short  periods  of  time — to  calculate 
out  sole  accurate  integral  cycle  precision  values.  Moreover,  the 
requirement  associated  with  continuous  measurement  sequences  of 
200s  or  more  is  that,  due  to  making  fast  approximations,  a  definite 
number  of  measurement  data  sequences  are  needed. 

The  fast  positioning  algorithms  introduced  by  this  article  do 
not  calculate  initial  integral  cycle  ambiguity  values.  Option  is 
made  for  the  use  of  measurement  sequence  window  type  computations. 
With  regard  to  each  observation  epoch,  use  is  made  of  forward 
shifts  associated  with  5  point  computation  windows.  Step  by  step, 
the  5  point  position  coordinates  associated  with  the  epoch  in 
question  are  calculated.  After  that,  the  average  value  for  the  5 
point  coordinate  solutions  acts  as  the  coordinate  for  that  epoch. 
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Fig.l  Observation  Epoch  5  Point  Position  Coordinate  Solution 

^  Solution  Window  (2)  Measurement  Sequence 

(3)  Measurement  Epoch 


The  calculation  process  is  as  follows. 

Take 

A<t>^  X(N[+ n[)+ <J>[- I 

to  act  as  carrier  wave  measurement  difference.  Take 

AA^'  = 


to  act  ^  as  initial  integral  cycle  ambiguity  difference .  In 

there  is  no  loss  of  lock  during  test  measurement 
processes,  the  values  in  question  do  not  vary  during  the  entire 
test  measurement  process.  With  regard  to  the  same  observational 
epoch,  equation  (6)  can  be  written  as  /24 

p;+  XAN‘+  A<J>=  [(r- (y-  y  V+  (z‘- AA/  I  (7) 


In  toe  equation,  unknown  numbers  include  ANi  ,  Xk  ,  Yk  ,  Zk  as  well 
as  Adri  and  ANi  and  are  reference  station  and  user  receiver  initial 
ambiguity  differences  with  regard  to  toe  ito 
satellite  and  take  the  form  of  whole  numbers. 

In  cases  where  there  is  no  loss  of  lock  during  test 
measurement  processes,  when  there  is  continuous  observation  of  a 
number  of  epochs,  toe  unknown  nximbers  ANi  do  not  vary.  However  Xk 
Yk  ,  Zk  as  well  as  Adri  have  some  differences  in  each  epoch. ' 
However,  with  regard  to  Adri  ^ 

Adr‘+  C(dT^-  dMJ+  (V^~  V^)  I 


In  tois  ej^ression,  as  far  as  C(drk  -  drr)  ,  that  is,  distance 
errors  given  rise  to  by  GPS  receiver  clock  errors  relative  to 
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reference  station  receivers  is  concerned,  this  term  can  also  be 
expressed  as 


C(Ar,+  AdrO^-  r„)) 


In  the  expression,  ATo  is  initial  observation  epoch  to  .  In 
respect  to  absolute  deviations  associated  with  instantaneous  GPS 
receivers  relative  to  reference  station  receiver  clocks,  during 
test  measurement  processes,  the  terms  in  question  are  maintained 
invariable.  Adt  is  the  degree  of  deviation  given  rise  to  by 
internal  clock  oscillator  frequency  instabilities  of  the  two 
receivers.  With  regard  to  the  quartz  oscillator  frequency  scale 
used  by  general  GPS  receivers,  the  stability  is  on  the  order  of 
magnitude  of  1010  .  For  this  reason,  Adt  should  be  higher  than  an 
order  of  magnitude  of  1010  .  Then,  the  changes  of  C(ATo  +  Adt(tto 
-  to  ) )  in  middle  and  low  dynamic  states  will,  in  all  cases,  not  be 
lower  than  the  order  of  magnitude  of  10-2  m/s.  Relative  to  LI  wave 
band  carrier  wave  wave  lengths  of  19.3cm,  it  is  possible  to  assiame 
that  there  is  basically  no  change  between  the  jth  observation  epoch 
and  the  j+lth  observation  epoch  Adt.  Thus,  it  is  possible  to 
assume  that,  as  far  as  C(dTk  -  drr)  associated  with  the  term  in 
question  is  concerned,  there  is  basically  no  change  between  the  jth 
observation  epoch  and  the  j+lth  observation  epoch. 

With  regard  to  (dMk  -  dMr) ,  due  to  the  fact  that  correlations 
user  GPS  receivers  relative  to  reference  station  receivers  are 
relatively  small ,  the  terms  in  question  are  important  sources  of 
errors  during  measurement  and  positioning.  As  far  as  concrete 
environments  with  different  positioning  landforms  are  concerned, 
there  are  also  some  differences  in  changes  between  two 
observational  epochs.  However,  with  regard  to  reference  station 
receiver  dMr  ,  there  is  basically  no  change.  Changes  between  two 
observational  epochs  are  fundamentally  changes  associated  with  user 
GPS  receivers  in  dynamic  environments.  Option  is  made  for  the  use 
of  specialized  antennas.  With  respect  to  static  states,  use  is 
made  of  (dMk  -  dMr)  between  two  observational  epochs,  and  there  is 
basically  no  change.  In  regard  to  dynamic  state  utilizations — in 
extreme  cases  the  terms  in  question  lead  to  errors  estimated  on 
the  order  of  magnitude  of  10-2m. 

With  regard  to  (Vk  -  Vr) ,  GPS  receiver  correlations  relative 
tc  reference  station  receivers  are  also  relatively  small.  However, 
the  quantities  in  question — in  certain  dynamic  utilizations — 
possess  a  definite  randomness.  Going  through  measurements 
associated  with  multiple  observation  epochs,  applications  of 
statistical  treatments  are  capable  of  the  influences  on  the  order 
of  magnitude  of  10-2m  . 

Thus,  as  far  as  Adri  between  each  observational  epoch  is 
concerned,  it  can  be  assumed  that  Adri  between  the  jth 
observational  epoch  and  the  j+lth  observational  epoch  is  basically 
invariable,  with  respect  to  deviations  between  Adri  channels,  they 
primarily  caused  due  to  noise  associated  with  various  receiver 
channels,  that  is,  they  are  still  under  the  influence  of  (Vk  -  Vr) . 
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sources  in  question  are  capable  of  being 
^  receiver  designs  and  statistical  processing  of  data 
obtained  in  measurments.  Pseudo  range  errors  given  rise  to  by  the 
S'^iSve??^  question-in  extreme  cases-are  capable  of  reachSg 

Co^ining  pseudo  range  changes  given  rise  to  by  all  the  error 
quantities,  the  error  grand  total  is  still  on  the  cm  level,  it  il 
possible  to  assume  that  Adri  between  various  channels  as  well  as 
obsOTational  epoch  is  invariable.  The  deviations 
manifested  in  differences  associated  with  initial  integral  cycle 
ambiguity  yalues  ANi  are,  thus,  selected  as  ^  cycle 


Adr  =  Adri 


/25 

US  making  obseryations  associated  with 
respect  to  4  satellites  is  concerned,  it  is 
possible  to  arriye  at  4  single  error  equations  of  the  type  of 

^  7k*  8  unknown  numbers: 

/  and  ANi  (i=l,2,3,4) .  There  is  no  way  to  solye 

flL  i?e  ^"succession,  2d  epoch  obseryations 

also  are  capable  of  arriving  at  4  single  error  equations  of  the 

o^ainfd^S^^Jh^  the  4  single  error  equations 

obtained  for  the  1st  epoch,  there  are,  then,  8  single  error 

e^ations.  However,  this  epoch  also  adds  3  unknown  numbers:  xk  , 

total  for  the  8  single  error  equations, 

Sdl?,™  following  along  with  the 

addition  of  observational  epochs,  each  epoch  is  capable  of  arrivina 
at  4  single  error  equations  of  the  type  of  equatiSi  (7)?  However 

sdded.  After  continuously  observing  5 
pochs,  it  IS  possible  to  reach  20  single  error  equations  of  the 

nSbers  time,  there  are  20  unknown 

Common  unknown  ntimbers :  Adr,  ANi  ( i=i ,2,3,4) 

Respective  epoch  unknown  numbers:  Xk  ,  Yk\  ^Zk  (k=l,2,3,4,5) 

It  is  then  possible  to  calculate  meaningful  solutions.  After 
continuing  with  observations  of  the  next  epoch,  they  are  then  set 

equations  of  Se  t?|e  of 

equation  p)  obtained  from  the  previous  4  epochs.  It  is  also 
possible  to  calculate  a  set  of  solutions.  Thus,  beginnina  with  5 
point  calculation  window  displacement  calculations — after  carrying 
out  5  Iterations  of  window  calculations— each  epoch  is  capable  of 
one  set  of  calculated  values 


X^Y^Z^kj  \ 

J  J  J  I 

1.  2.  3,  4.  5)i 
0=  1.  2.  3.  4.  5)  * 
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*  y-  *  yk 

Taking  the  5  coordinate  values  ^  ’  j  '  ^  j  which  are 

obtained  in  association  with  calculations  for  each  observation 
epoch,  option  is  made  for  the  use  of  solutions  for  mean  values  to 
obtain  positioning  solutions  Xk  ,  Yk  ,  Zk  . 

4  CONCLUDING  REMARKS 

Through  analysis  with  regard  to  Adr — as  far  as  option  being 
made  for  the  use  of  carrier  wave  phase  difference  positioning  is 
concerned — residual  errors  are  already  very  small.  With  regard  to 
positioning  solutions  obtained  with  the  calculation  methods  in 
question — due  to  opting  for  the  use  of  multiple  point  average  wave 
filtering~to  a  certain  degree,  it  is  possible  to  reduce  errors 
given  rise  to  by  receiver  internal  noise.  In  the  end,  positioning 
precisions  are  able  to  achieve  centimeter  levels. 
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